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We analyze a traffic overflow system that consists of two groups of 
trunks, with waiting spaces for each group, and some overflow ca- 
pability from the primary to the secondary group. We consider the 
case in which the number of waiting spaces in the primary queue is 
large compared to the corresponding number in the secondary queue 
and to the number of trunks in the secondary group. The case of an 
infinite number of waiting spaces in the primary queue is also 
allowed. We contrast the approach presented with some previous 
approaches that are suitable when the number of waiting spaces in 
the primary queue is not comparatively large. As with previous 
approaches, the aim is to reduce the dimensions of the system of 
equations to be solved in order to calculate various steady -state 
quantities of interest. Our results include expressions for the loss 
probabilities, the probability of overflow from the primary to the 
secondary group, and the average waiting times in the queues. We 
also obtain the stability condition under which the results are valid 
when the number of waiting spaces in the primary queue is infinite. 

I. INTRODUCTION 

In this paper, a particular traffic overflow system with queueing is 
analyzed. The same system has been investigated previously where 
techniques were developed for reducing the dimensions of the system 
of equations to be solved in order to calculate various steady-state 
quantities of interest. 1,2 This analysis is a considerable improvement 
over the earlier study for the case when the number of waiting spaces 
in the primary queue is large compared to the corresponding number 
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in the secondary queue and to the number of secondary trunks. 1 A 
suitable analysis has already been made for the case when the number 
of waiting spaces in the secondary queue is large compared to the 
corresponding number in the primary queue and to the number of 
trunks in the two groups. 2 

In cases of interest, the sparse system of linear equations for deter- 
mining the steady-state probabilities of the number of demands in the 
two groups may have very large dimensions. The reduced system of 
equations then has considerably fewer dimensions, but the system will 
be dense. Numerical results based on the solution of the original sparse 
system of equations by successive over-relaxation techniques, and on 
the solution of the reduced system of equations obtained earlier by 
Morrison, 1 were presented by Kaufman, Seery, and Morrison. 3 The 
numerical values of the various steady-state quantities of interest 
obtained by the two procedures agree to many significant figures. 
Other procedures for solving the original system of equations are 
discussed in Refs. 1 and 4. 

The dimensions of the reduced system of equations obtained in this 
paper are independent of the number of waiting spaces in the primary 
queue. On the other hand, the dimensions of the reduced system of 
equations derived in Ref. 1 continue to increase as the number of 
waiting spaces in the primary queue increases. This is even more 
drastically the case for the original sparse system of equations. Con- 
sequently, the reduced system of equations that we derive is advan- 
tageous when the number of waiting spaces in the primary queue is 
large, and even more so when the number is infinite. 

The traffic overflow system considered consists of two groups, a 
primary and a secondary, with n k > servers and q k waiting spaces, 
which receive demands from independent Poisson sources S* with 
arrival rates \ k > 0, k = 1 and 2, respectively, as shown in Fig. 1. The 
service times of the demands are independent and exponentially dis- 
tributed with mean service rate /x > 0. If all n 2 servers in the secondary 
group are busy when a demand from S 2 arrives, the demand is queued 
if one of the q 2 waiting spaces is available; otherwise it is lost, that is, 
blocked and cleared from the system. Demands waiting in the second- 
ary queue enter service, in some prescribed order, as servers in the 
secondary group become free. 

If all /ii servers in the primary group are busy when a demand from 
Si arrives, and there is a free server in the secondary group and no 
demands waiting in the secondary queue, the demand is served in the 
secondary. If there are no free servers, then the demand is queued in 
the primary queue, if one of the qi waiting spaces is available; otherwise 
it is lost. Previously, two different cases were considered for the 
treatment of demands waiting in the primary queue. 1,2 In case I, a 
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Fig. 1 — Mean flow rates for an overflow system with queueing. 

demand waiting in the primary queue may enter service either in the 
primary group, when a server becomes free, or in the secondary group, 
if a server becomes free and there are no demands waiting in the 
secondary queue. In case II, no overflow is permitted from the primary 
queue, so that a demand in the primary queue must wait for a server 
in the primary group to become free. 

We consider only case I and use a different approach to analyze the 
overflow system. The new approach is preferable to the earlier ap- 
proach 1 when <7i is large compared to q-i and n 2 , or even infinite. The 
alternate earlier approach 2 is preferable when q* is large compared to 
<7i, Tii, and ri2, or even infinite. 

Let pu denote the steady-state probability that there are i demands 
in the primary group and j demands in the secondary group, either in 
service or waiting. These probabilities satisfy a set of generalized birth- 
and-death equations, which take the form of partial difference equa- 
tions connecting nearest neighboring states. The basic technique is to 
separate variables in regions away from certain boundaries of the state 
space, the elements of which are (i,j). These regions are shown in Fig. 
2a. The analogous regions corresponding to the earlier analyses 1 ' 2 are 
shown in Figs. 2b and 2c. The separation of variables leads to two sets 
of eigenvalue problems for the separation constant. The eigenvalues 
are roots of polynomial equations. The probabilities py are then 
represented in terms of the corresponding eigenfunctions. 
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The constant coefficients in these representations are determined 
from the boundary conditions and the normalization condition that 
the sum of the probabilities is unity. There are two sets of constants, 
corresponding to the representations of the probabilities py in the two 
shaded regions in Fig. 2a. In general, these constants have to be 
evaluated numerically. However, we do show that the set of constants 
corresponding to the representation in the primary queueing region 
may be expressed explicitly in terms of the other set of constants. This 
further reduction leads to a problem of fewer dimensions than the one 
obtained in the first analysis, 1 regardless of the size of q\ > 0. 

Various steady-state quantities of interest may be expressed in terms 
of the probabilities p./. Among these are the loss (or blocking) proba- 
bilities, the probability of overflow from the primary to the secondary 
group, the probabilities that a demand is queued, and the average 
waiting times in the queues. These quantities may be expressed directly 
in terms of the constant coefficients occurring in the representations 
of the probabilities pi y , averting the need to calculate the pij. This is a 
key advantage of the reduction in the dimensions of the problem. 

The solution is built up in stages. In Section II, we consider solutions 
to the birth-and-death equations in the region < i < zii — 1, < / < 
/i2 + q 2 of state space. This corresponds to the case q\ = 0, and was 
analyzed earlier as a special case. 1 Its inclusion here is for the sake of 
completeness. In Section III, we examine the solution in the region of 
state space corresponding to queueing in the primary, and give a 
heuristic derivation of the stability condition when q\ = oo. Section IV 
discusses the boundary and normalization conditions, while Section V 
is devoted to the calculation of various steady-state quantities of 
interest. In Section VI, we show how to achieve a further reduction in 
the dimensions of the problem by the introduction of a generating 
function. This reduction is obtained in another way in Section VII. 
Properties of the eigenfunctions that occur in the representations of 
the probabilities p iy are given in Appendix A, the eigenvalues corre- 
sponding to the primary queueing region are discussed in Appendix B, 
and results pertaining to the generating function are derived in Appen- 
dix C. 

II. REPRESENTATION OUTSIDE THE PRIMARY QUEUEING REGION 

For convenience, we define 

ki = ni + q u k 2 = n 2 + q 2 . (1) 

Let pij (0 < i < k\, <j < k 2 ) denote the steady-state probability that 
there are i demands in the primary, and j demands in the secondary. 
Define also the traffic intensities 

en = Xi/n, a 2 = A 2 //t. (2) 
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As usual, we let 6y denote the Kronecker delta, i.e., 

*-{j £{ <3) 

Finally, we introduce the function 

_fl» ^o, (4) 

*"{<>, Z<0. (4) 

The probabilities py satisfy a set of generalized birth-and-death 
equations. For case I, it may be shown 1 that 

[aid - Sik.xj-^) + a 2 (l - 8jk 2 ) + min(i, /ii) + min(y', n 2 )]pij 

= (1 - X'-»i-iX"2-y-i)[ai(l ~ 5 «o)p«-i,y 

+ (1 - 8,* 2 )minO' + 1, n 2 )pij + i\ (5) 

+ (1 - MfaiSm.X^-T + a 2 (l - X«-n,-lXn2->)]P'.>-l 

+ (1 - 6 flkl )[(l - x.-n,Xn 2 -y-i)niin(i + 1, rn) + n 2 Xi-n&»*]Pi+\,j, 

for < i < ki, <7 < £ 2 . These equations were constructed to imply 

that 

Pij = 0, /ii + 1 < i < Aii, <y < n 2 - 1, (6) 

since there can be no queued demands in the primary when there is a 
free secondary server. 
The normalization condition is 

,=0 y-o »-n,+l y'-n^ 

For < i < /ii - 1 and < /' < & 2 , the variables in (5) may be 
separated, 1 and there are solutions of the form a,/?,, where 

{ax + i + p)a, = ai(l - 5«o)a«-i + (*' + Dam, (8) 

for < i < n\ — 1, and 

[a 2 (l - 6/0 + min 0> "2) - p]pj 

= a 2 (l - 8j )Pj-i + (1 - 5y* 2 )min(y + 1, re 2 )j8/+i, (9) 

for < 7" < A> 2 , and p is a separation constant. The solution of (8) 
may be expressed in terms of Poisson-Charlier polynomials. 5 ' 6 We 
denote here the solution of (8) for which a = 1 by Si(p, Oi). The prop- 
erties of Si{\, a) which we will need in the analysis are contained in 
Appendix A. 
For </ < 7i 2 - 1, (9) becomes 

(a 2 +/ - p)fy = 02(1 - S j0 )fij-i + (j + l)$j+i, (10) 
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and we see that there are solutions fij proportional to Sj(—p, a 2 ) for 
<7 < n 2 . For n 2 <y < k 2 , we find that 



[a 2 (l - 8jk 2 ) + n 2 - p]fij = a 2 pj-i + n 2 (l - Sjkjfy+i. 
To find a representation for solutions to (11), we define 1 



MM,.l-|j|| 



//2 



c/i 



fite + 712 — p 



02/12 



and 



hip) = *k 2 -j(p) - Vkt-j-M, 



(11) 



(12) 



(13) 



where Ui(x) denotes the Zth Chebyshev polynomial of the second 
kind. 7 The properties of these functions are discussed in Appendix A. 
However, we note here that Uo(x) = 1, U-i(x) = 0, and hence 
<f>* 2 (p) ■ 1. From (11), and (134), it follows that /?> is proportional to 
hip) for/i2 - 1 <j<k 2 . 
Thus, as a basic solution to (9), we take 



fij 



sj(~p, a 2 )<j>n 2 (p), 
Sn.X-p, a 2 )<t>j(p), 



0<y'</i2, 

n 2 - 1 <y'< k 2 , 



where 



Sn 2 -i(-p, a 2 )<j> n2 (p) = s„^(-p, a 2 )<>^-i(p). 



(14) 



(15) 



With the help of (13), (128), (129), and (133), (15) may be expressed in 
the form, 

p[ Sn2 (l - p, a 2 )¥ 92 (p) - Sn 2 -i(l - p, a^^-xip)] = 0. (16) 

The expression in the square brackets is a polynomial in p of degree 
k 2 . It was shown 1 that its zeros are positive and distinct, and we denote 
them by p m , m = 1, 2, • • • , k 2 . We adopt also the convention po = 0. 
The calculation of the {p m } for m > is related to an eigenvalue 
problem for a symmetric tridiagonal matrix. 1 

For the moment, we consider q\ = 0, so that k\=n.\. This corresponds 
to the case of queueing in the secondary only, which has been analyzed 
as a special case. 1 In terms of the separated solutions to (5) for 
< i < n\ — 1, a representation of the solution for < i < n\ , <j < k 2 
is 



Po = 



2 d m Si(p m , ai)Sj(-p m , a 2 )<t>n 2 (pm), < j < 7l 2 , 



£ d m Si(p m , a^s^i-pm, a 2 )<t)j(p m ), n 2 £7 < k 2 . 

m— 



(17) 
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The constants {d m } are determined only to within a multiplicative 
factor by the boundary conditions at i = n h < / < k 2 , since the 
boundary condition at / = n 2 may be shown to be redundant. 1 The 
multiplicative factor is determined by the normalization condition (7). 
In general, the constants {d m } have to be calculated numerically. 

For q x > 0, the representation (17) may still be used for < i < n it 
< j < k 2 , but a different representation must be found for /ii + 1 < 
i < k u n 2 < / < k 2 . The earlier approach 1 was to find separated 
solutions to (5) in the region n x + 1 < i < k x , n 2 + 1 </ ^ k 2 , i.e., away 
from the boundary at / - n 2 . This led to a representation of />,> for 
ni<i<ki,n 2 <j<k 2 . Solutions were found that match those in (17) 
for i = ni t n 2 <j <k 2 , and an additional q\ solutions were found that 
vanish at i = /ii . By the superposition of a suitable linear combination 
of these q\ solutions, the boundary conditions at/ = n 2 , n x + 1 < i < 
k\ could be satisfied, without disturbing the matching of the represen- 
tations at i = n u n 2 <j<k 2 . The constants {d m } are still available to 
satisfy the boundary conditions at i = m, < / < k 2 , and again the 
condition at/ = n 2 is redundant. 

Since the number of additional constants to be determined in the 
above approach is q x , this procedure is not suitable if q\ is large. In the 
next section, we investigate an alternate approach that is suitable if q x 
is large or even infinite. 

III. REPRESENTATION IN THE PRIMARY QUEUEING REGION 

We now assume that qi > 2 and q 2 > 1, although these restrictions 
may be relaxed. For n a + 1 < i < k x - 1 and n 2 <j<k 2 , the variables 
in (5) may be separated and there are solutions to the partial difference 
equations of the form 



% 7li<l'<A5l, n 2 </</2 2 , 



where 



n x 



a 



- — (1 - a) + a 2 (l - 8jk 2 ) + n 2 (l - 8 jn2 a) 



fj 



= a 2 {\- Bjnjfj-i + n 2 (l - Sfljfm* n 2 <j<k 2 , (18) 



and a is a separation constant. Since cti ^ 0, and we are interested in 
nontrivial solutions, a ^ 0. In particular, for n 2 + 1 5/ ^ k 2 , we have 

[a 2 (l - «>*,) + n 2 - w]/; = a^-i + n 2 (l - Sjijfj+i, (19) 

where 

(20) 



w = I — — Hi 1(1- a). 
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We see from (134) that a solution which satisfies the boundary condi- 
tion at j = k 2 is 

fj = f k2 <t>Au), n 2 <j<k 2 , (21) 

with <£, defined as in (13). 
It remains to satisfy the boundary condition aty* = n 2 , namely, 

[a 2 + n 2 (l - a) - w]/^ = n 2 fn 2 +i. (22) 

With the help of (21) and (13), we obtain 

mil - a)4h^(a) + (a 2 - «)[¥, 2 («) - % 2 -du>)] 

- /12^-iM - % 2 -2(u)] = 0. (23) 
This reduces with the help of (20) and (133) to 

(1 - a)H(a) m 0, (24) 

where 

H(a) = lrn- ^W 92 (w) + n^w), (25) 

with (o defined as in (20). Hence, a = 1 or H(a) = 0. 

Now ^/(co) is a polynomial of degree I in u. Hence, a q2+1 H{a) is a 
polynomial of degree 2q 2 + 1 in o, so that H(a) = has 2q 2 + 1 roots. 
In Appendix B, we show that these roots are positive and distinct, and 
that at least q 2 of them lie between and 1, and at least q 2 of them are 
greater than 1. More precisely, define 



A\ = n.\ + n 2 



* g2 -i(0) 
** 2 (0) 



(26) 



Then, if cti <A U q 2 + 1 roots lie between and 1, and q 2 roots exceed 
1. If a.\ > A\t then q 2 roots lie between and 1, and q 2 + 1 roots exceed 
1. Further, if ai = A u then H{\) = 0. We denote the roots of H(a) = 
by a r , r = 1, • • • , 2q 2 + 1, withO < o\ < — < o 2q2 +\ . For convenience, 
we define ao = 1. Then a Q2 < 1 and a q2+2 > 1, and a 92+ i < 1 if ai < Ai, 
and <jq 2 +i > 1 if a\ > A\. Note that a r , r = 0, • • • , 2q 2 + 1 are distinct 
if d 5^ Ai . On the other hand, ao = a q2 +\ if a\ = A\. 

For q\ < oo, we assume for simplicity that a\ ^ Ai. We then have 
2g2 + 2 distinct roots, and it follows that we may represent the 
probabilities in the form 

2 92 +l 

Pij= £ e^ r ~ n '<$>jM, /I, <i<Ai, n 2 <j<k 2 , (27) 
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where 



<O r = I Hi 1(1 - Or) 

and the constants {e r } are to be determined. Consistency of the 
representations (27) and (17) at i = m requires that 

2q 2 +l *2 

J er^y(tOr) = £ d m S n ,(pm, Ct^S^-pm, a 2 )<t>j(Pm), (28) 

r-0 m-0 

for n 2 <j<: k 2 . 

For gi = oo, we assume that Oi < Aij the reason for this will be 
apparent shortly. We use the representation (27), but the normalization 
condition (7) implies that e = and e r = 0, r = q 2 + 2, • • • , 2g 2 + 1 
since a = 1 and a r > 1, r = q 2 + 2, ■ . . , 2q 2 + 1. We thus have q 2 + 1 
conditions to replace the q 2 + 1 boundary conditions lost on removal 
of the boundary at i = ki. If a\ > A u then a 92+ i > 1, and we would 
have to take e q2 +i = also. However, this would leave only q 2 constants 
e r , r = 1, • • • , q 2 to satisfy the q 2 + 1 conditions (28), which is 
insufficient. 

The stability condition a a < Ai, for gi = oo, may be derived heurist- 
ically as follows. Consider the situation when there are demands 
waiting in the primary queue. It follows that all n 2 servers in the 
secondary are busy. As the point of instability is reached, the primary 
queue will never be empty, and hence acts as an infinite source for the 
secondary servers. Demands from this source enter the secondary only 
when a server becomes available and there are no demands waiting in 
the secondary queue. Since the secondary system is a loss system, it is 
ergodic. 

Let Pj,j = n 2 , • • • , k 2 , be the steady-state probability that there are 
j demands in the secondary. Then the rate at which demands in the 
primary queue are served in the secondary is mfiP^, and the rate at 
which they are served in the primary is #ii/i. Since Xi = ai/i f it follows 
that the stability condition is 

ai < /ii + 7t2iV (29) 

But, as for the system without the infinite source, 8 

a 

P^, n 2 ^j<k 2 . (30) 



■<-© 



With the infinite source present, P, = 0, > = 0, • • • , n 2 - 1, and P^ is 
determined by 

£P, = 1. (3D 
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From (30), (31), and (136), we obtain 

*, 2 -i(0) 



Pn,= 



1 - 



*„(0) 



(32) 



From (26), (29), and (32), it follows that the stability condition for 
q\ = oo is ai < A\. 

IV. BOUNDARY CONDITIONS 

The constants {d m } and {e r } introduced in (17) and (27) must be 
chosen so as to satisfy the boundary conditions at i = /ii and i = k\ . 
The boundary conditions at i ■» k\ imply that 

[a 2 (l - Sjk 2 ) + m + n 2 ]pk i j 

= a\pk x -\j + 02(1 — Bjn t )pk l ,j-\ + n 2 {\ — Sjkjp^j+i (33) 
for n 2 Sj < k 2 . At i = »i, we have 
(ai + a 2 + m +j)p ni j 

= aip„,-i.> + (ai + a 2 )(l - S/o)P»,.y-i + + l)p»,./+i (34) 
for0<y* <n 2 - 1, 
(d + a 2 + /ii 4- n 2 )/) n ,^ 2 

= aiPm-1^2 + (Ol + a 2 )Pn,.n 2 -l + »2/>n,^2+l + (ll + /l2)Pn,+1^2, (35) 

and 

[ai + a 2 (l - 5/* 2 ) + n\ + n 2 ]pn vJ 

= aiPn,-l.> + a 2 p ni J-l + /llPn,+l,y + W^l ™ Sjh^Pn^J+l, (36) 

for n 2 + 1 <y < £ 2 . 

If we substitute (17) in (34), we find, on reduction with the help of 
(125), (128), and (129), that 

2 dm<f>n 2 (pm)[PmS„,(l + P«i «l)*>( — P«t ^ 2 ) 
m-0 

+ aiSi,,(pm, ai)s,(-l - p m , a 2 )] = 0, (37) 

for <y < /i 2 - 1. Similarly, if we substitute from (17) and (27), (36) 
implies, after reduction, that 

£ d m p m s„,(l + p m , ai)s„ 2 (-p m , a 2 )<f>,(p m ) 

m-l 

2*72+1 

- £ e r (nia r - ai)<M« r ) = 0, (38) 

r-0 
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for n 2 + l<j<k 2 , and (33) implies, with the help of the fact that fj = 
hM satisfies (18) for < r < 2q 2 + 1, that 

2q 2 +l 

£ e r oV[(ni + n 2 8 jn2 )a r - a,]</> y (w r ) = 0, (39) 

r-0 

for n 2 <y < k 2 . Finally, (17), (27), (15), (125), (128), and (129) in (35) 
imply that 

A 2 

£ d m <S>n 2 {Pm)[pmSn^ + Pm, Ct^S^-Pm, d 2 ) 

m-Q 

+ ais„,(p m , ai)s„ 2 (-l - p m , 0,2)] 

2q 2 +l 

-(/ii + n 2 ) 2 erO^M = 0. (40) 

r-0 

We will now show that the boundary condition at i = /ii , j = n 2 is 
redundant and may hence be dropped. Summing (39) from j = n 2 to 
k 2 , we find that 

2<7 2 +l 

£ e r a?'[(nia r - ai)^ 92 (w r ) + n 2 a^(co r )] = 0. (41) 

r-0 

For r = 1, 2, . • • , 2q 2 + 1, it is true that 

if( CTr ) ■ I m - — \* Q2 M + nrf^M = 0, (42) 

and ffo = 1 and <o - 0, so that (41) and (26) imply that 

e Q (Ai - a\) = 0, 

that is, 

e = 0, (43) 

since ai ^ A\ by assumption. The redundancy of the boundary condi- 
tion (40) may then be verified by summing (37) fromy" = to n 2 - 1, 
(38) from./' = n 2 + 1 to k 2 , and (28) from; = n 2 to k 2 , and by using the 
facts that p m satisfies (16) for m = 0, 1, 2, • • • , k 2 , and e = 0, together 
with the recurrence relations found in Appendix A. 

Consider (39) for./ = n 2 . We note that <>„ 2 (w r ) is nonzero; otherwise 
(22) with fj = <J>,(co r ) would imply that<f> n2+ i(w r ) = 0, and consequently 
that tt>A<0r) = 0, in contradiction to the fact that <f>k 2 M = 1. Further, 
since the roots of if (a) are distinct, it follows that not all of the square- 
bracketed terms are zero. Hence, the equation is nonvacuous. Since 
summation of (39) fromy" = n 2 to k 2 implies e = 0, if we set e = in 
(38) and (39), then (39) fory* = n 2 is redundant and accordingly may be 
dropped. We remark that even for a,\ = Ai, we are still free to take 
e = since a g . i+ i = 1 in that case. As before, the boundary conditions 
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at y = n 2 (i = ti\ and i = &i) are redundant. Therefore, it follows that 
the constants {d m } and {e r } are determined by the boundary conditions 
at i = ri\ for < j < n 2 — 1 and n 2 + 1 < J : < k 2 , the consistency 
relationships at i = /ti, n 2 <y < A 2 , the boundary conditions at i = ki, 
n 2 + 1 </' < £ 2 if ^i is finite, and the normalization condition (7). 
We summarize the results of the last three sections as follows: 



Pij = 



m-0 

for < i < ni and, 



*2 

J d m Si(p m , (li)Sj(-p m , a 2 )<t>n 2 (pm), <> < rt 2 , 

m-0 

j d m Si(p m , ai)s^(-p m , a 2 )#/(p m ), /l2 <y* < A 2 , 



(44) 



2^2+1 



Pv - £ e^ "'^(air), n 2 <j<k 2 , 



(45) 



for Tii < i < Ai, (<7i < oo), where 

Pm[s^(l ~ pm, a2)^q 2 (pm) ~ S^_i(l - p m , 2 )^, 2 -i(p m )] = 0, (46) 

for m = 0, 1, 2, • • • , k 2 and 

tf(a r ) =0, r = 1, 2, 3, • ■ • , 2q 2 + 1 



(47) 



ai 



With 03 r = I All | (1 — CTr). 

The ^2 + 2q 2 + 2 constants {d m } and {e r } are determined, with eo = 0, 
by (28) and (37) to (39) only to within a multiplicative constant, which 
is determined by the normalization condition. From (44), (46), (13), 
(130), it follows that 

*2 

£ Pij = d s,(0, ai)[sn 2 d, a 2 )* Q2 (0) - s^-id, a 2 )¥ 92 -i(0)], (48) 
y-o 

for < i < ni, so that (7), (13), (45), and (130) imply that 
dos„,(l, ai)[s„ 2 (l, a 2 )^ 2 (0) - s„ 2 -,(l, a 2 )¥ 92 _i(0)] 

292+1 



+ X e r r 9l (o r )^ (72 ( Wr ) = 1, (49) 

r=l 



where 



iya - V r -1 = 



,-/ _ j *a - r>/u - a, s*i, 



</, 



1 = 1. 



(50) 



We point out that the special case gi > 1, ^ 2 = has previously been 
analyzed. 1 In the primary queue region, solutions of the form 



Pi,n 2 = 



a\ 



n,i + n 2 



p nx ,n 2 , nx<i<ku (51) 
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where 

P»,.ii, = S d m s ni (Pm, a^s^i-pm, a 2 )<K(P"»)» (52) 



m-0 



were found. We note that since k 2 = «2 for q 2 = 0, <f>« 2 (p) = <f>k 2 (p) ■ 1, 
above. For <ji = », the stability condition was found to be a\ < rt\ + n 2 . 
It is readily verified that the representations (44) and (45), and the 
stability condition d < A\, produce the same result. 

For <7i = oo and q 2 5: 1, we must set e r = for r = q 2 + 2, • • • , 
2q 2 + 1, as well as e = 0, as discussed in Section III. The remaining 
k 2 + q 2 + 2 constants {d m } and {e r } are then determined by (28), (37), 
and (38) to within a multiplicative constant, which is determined by 
the normalization condition. But, from (50), 

r-(0 = Hm T qi (& = TT-^> < Z < 1. (53) 

Hence, for q\ = oo, the normalization condition (49) is 
cfo8„,(l, Oi)K(l, a 2 )¥ 92 (0) - Sn 2 -ld, O2)*„-l(0)] 



q 2 +l 



e r o r 



+ 2 TT^-T^W-l. (54) 

r _0 (1 — Or) 

V. SOME STEADY-STATE QUANTITIES 

We now proceed with the calculation of various steady-state quan- 
tities of interest. These quantities are shown in Fig. 1, which depicts 
mean flow rates. The loss probabilities are, for the primary and 
secondary, respectively, 

Lx - 2 pjw. L * - 2 Pi** < 55 > 

while the probabihties that a demand from the primary, or secondary, 
source is queued on arrival are 

Qi = Y 2 J*. (9i 2= 1), Q 2 =H Pifc (92 ^ 1). (56) 

i-n, y-nj *'-0 /-"| 

The probability that a demand arriving from the primary source 
overflows immediately is 

hi = "j Pn,./. (57) 

y-0 

Since the mean service rate is fi for each server, the mean departure 
rate from the primary queue to the primary servers is 

flu = #n/i 2 2 PV> (9i S D, ( 58 ) 

i-n,+l y'-n^ 
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while the mean departure rate from the secondary queue is 

#22 = n 2i u £ I Pv, (?2 > 1). (59) 

i-0 7-/I2+I 

The mean rate at which queued demands in the primary overflow to 
the secondary servers is 

*. 
Ri2 = n 2 fi £ p,,„ 2 , (?i>l). (60) 

i-n,+l 

The average queue populations are, for the primary and secondary, 
respectively, 

V,- t £ (i-ni)pv, (91 >1), (61) 

and 

V 2 =£ £ U-n 2 )pij, (q 2 >1). (62) 

Also, the average number of demands in service in the two groups are 
Xx = £ £ ipij + Hid " Vo) £ £ PU> (63) 

i—O 7— i— t^+1 7— /lj 

and 

X2 = £ £ jjty + n 2 £ £ Pv- (64) 

1=0 j— 1=0 y— nj 

An application of Little's theorem 8 to the primary and secondary 
queues shows that the average waiting times of the queued demands 
in the primary and secondary are given by 

Wi-^,(qi*D, W 2 = ^r-,(Q2^D, (65) 

respectively, independently of the service order within each queue. 
Also, if we apply Little's theorem to the primary and secondary groups 
of servers, we obtain 

Xi(l-Li-Ii 2 )-Ri2 = iiXi, A 2 (l - L 2 ) + Ai/ 12 + Ri 2 = nX 2 , (66) 

since the mean service rate is /x. 

The steady-state quantities of interest may be expressed directly in 
terms of the constants {d m } and {e r }, with the help of the representa- 
tions (44) and (45). From (55), with the help of (13), we find that 

2q 2 +l 

Li= £ e r o?* g2 M, (67) 
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and with the help of (130) and the fact that<fr 2 (p) ■ 1, we find that 

k 2 292+1 

L 2 = £ d m Sn x {\ + Pm, a^Sn^-Pm, a 2 ) + £ e r r Ql {a r ), (68) 

m-0 r=1 

where T q (£) is as defined in (50). Similarly, it follows from (57) and 
(130) that 

*2 
il2 = E dmSn^Pm, ai)8n,-l(l ~ pm, 02)<>n 2 (pm). (69) 

m-0 

Next, we find from (56), (50), and (13) that 

292+1 I p \ 

Qi = I (f J r g ,(a r )^ 2 (co r ). (70) 

Also, (58), (50), and (13) imply that 

2 92 +l 

i*n = /ii/x X dVa^W, (71) 

while from (60) and (50) it follows that 

2?2+l 

Ri2 = n 2t i t toTtiMfrJM- < 72 ) 

r=l 

It may be shown directly from (5) that 

fin + R12 = \iQu (73) 

and, using (42), it is readily verified that this is consistent with the 
representations (70) to (72). 
At this point, we define 

rj=Yp iJ , n 2 <y<A!2. (74) 

i=0 

If we sum on i in (5), we obtain 

[a 2 (l - 8jk 2 ) + n 2 ]rj = a 2 rj-i + n 2 (l - 8jk 2 )r J+1 , (75) 

for n 2 + 1 <j : < k 2 . It follows that 

rwj = a 2 rj- h n 2 + l<j <k 2 . (76) 

Hence, since L 2 = r k2 , from (55) and (74), we find that 

r,= (-W n 2 <j<k 2 . (77) 



a 2 
Then, it follows from (50) and (59) that 



fi 2 2 = A 2 rJ-|L 2 . (78) 



t22 — A 2 1 q 2 

Ml 
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Also, from (56) and (77), we deduce that 

Q 2 = r *(ijW (79) 

These representations satisfy R 22 = A2Q2. We remark that this latter 
result, as well as (73), hold since, in the steady state, the average 
arrival and departure rates are equal, for both the primary and sec- 
ondary queues. 
From (61) and (13), we find that 

2<7 2 +l 

Vi= X e r A Qx (or)% 2 M> (80) 

r-\ 

where 

A ^ fi* 1^(7 + 1), *«1. (81) 

We note in passing that 

r,tf ) - g^ 1 + (1 - £)A 9 (£), (82) 

a formula that may be numerically useful. From (62), (74), (77), and 
(81), it follows that 

Finally, from (63), with the help of (7), (48), and (131), we find that 

X, = m - dos„,-i(2, Oi)[«^(l, a 2 )* 92 (0) - s^-id, a 2 )¥ 92 -i(0)], (84) 
and, from (64), with the help of (7), (17), (130), and (131), 

X 2 = 712 - £ dmS ni (l + Pm, ai)Sn 2 -l(2 ~ p m , a 2 )<l>n 2 (Pm). (85) 

m-0 

We remark that in view of (67) to (69), (72), (84), and (85), the relations 
presented in (66) provide a useful numerical check. 

We point out that the results of this section up to this point are 
valid for qi < oo. We recall that for q\ = oo, e r = for r = q 2 + 2, • • • , 
2q 2 + 1, since a r > 1 for these values of r. In addition, the stability 
condition a\ < A\, with A\ defined in (26), is necessary. Under these 
assumptions, since < a r < 1 for r = 1, ■■■ ,q 2 + l,we find, analogously 
to (67), that 

LiL-. = 0, (86) 

as expected. Also, with the help of (53), we obtain, analogously to (68) 
and (70) to (72), 
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92 +l 



e&r 



L 2 L-.= 2 rf m Sn,(l + Pm, ai)Sn 2 (-pm, a 2 ) + 2 n " " y \ 81 > 

m-0 r-l U Or) 

and 



92+1 



e r a r 



«n|, 1 --ni/i 2 77^ **<«*), (89) 

r-l U ~~ a rJ 



and 



92+1 p /j 

r-l U "~ °r) 

Similarly, from (81), we have 

A„(0 = Urn Afctt) = -^4^ ' < « < L (91) 

Hence, corresponding to (80), we obtain 



92+1 



e r o r 



^"IFTF*^ (92) 

The expressions for Ji 2 , R 22 , Q2, Vt, X u and X 2 , given by (69), (78), 
(79), and (83) to (85), still hold. Of course, the constants {d m ) and 
{e m }, which occur in these expressions, (87) to (90) and (92), are those 
pertaining to q\ = 00, and are determined in the manner discussed at 
the end of Section IV. 

VI. ELIMINATION OF ONE SET OF CONSTANTS 

We will now show how to determine the constants {e r } in the 
representation (27) of the probabilities py for n x < i < k u n 2 <j < k 2 , 
in terms of the constants {d m ) in the representation (17) of the 
probabilities p tj for < i < Hi, <j < k 2 . For m < i < k u n 2 <> < k 2 
we have, from (5), 

[ai + a 2 (l - SjkJ + Mi + n 2 ]pij 

= aipi-u + [(ai + a 2 )8in 1 8jn i + a 2 (l — QjhJIPu-* 

+ (iii + n 2 8 jn2 )pi+ij + n 2 (l - Sjhjpij+i. (93) 

We consider this equation for m < i < °°, and introduce the generating 
function 

Gj(z) = 2 J** 1 -- (94) 

The series will be convergent for sufficiently small \z\. 
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We consider z ^ and define 



r-<i-*)(^-«i). 



(95) 



Then, from (93) to (95), it follows that 



a-Al-SjkJ l * 2 (l-j**) -{ 



G,(z) 



- a 2 (l - ty*)G/-i(ff) - n 2 (l - WW*) 

= aiPn,-l,y + (Ol + a 2 )8 jn2 Pn,J-l (»1 + I^Jn^P^J, (96) 

z 

for n 2 <j<k 2 . We now define Gj'^z) for n 2 < Z < A 2 as the solutions 
of the equations 



cf<*) 



Ob(1 - *>*,) + ns !--«>», )-r 



- 02(1 - W«i-i(«) - "2d - SjkJGfUz) Sjt, (97) 

for n 2 <y < A 2 . Then, from (96) and (97), we have 

*2 

G y (z) = X (nipn,,/ - a lZ p n ,-u)Gf\z) 

+ [/i2p n ,^ 2 - (ai + a 2 )2p„,^_i]GJ n2) (2). (98) 



We define 



>^2 



a 2 z 



(99) 



where ¥/(p) is as defined in (12). It is shown in Appendix C that the 
solution of (97) is 

IT, , I ( I I I 7 1 

n 2 <>< /, 



id) 



orw 




*i(OTy(«) 



(l-z)H(l/z)' 



/M*""* T,(z)fc(tt 



(100) 



(l-z)H(l/z)' 



l<j<k 2 , 



for n 2 < / < £ 2 , where <f>,(p) and if(a) are as defined in (13) and (25), 
respectively, with co as defined in (20). We know that H(a r ) = 0, r — 
1, • • • , 2<7 2 + 1, with < o\ < • • • < a 2<72 +i . As discussed in Section III, 
Oq 2 < 1, Oq 2+2 > 1, and a g2+i < 1 if d\ < Ai, and a, 2+ i > 1 if d > At, 
where Ai is given by (26). We assume that an*Ai t so that a Q2 +i ^ 1. 
Now, from (12), (13), (20), (25), (95), and (99), z q *H{l/z), *Mfo(fl 
and z'^Tjiz) are polynomials in 2 of degrees 2g 2 + 1. 2(k 2 — j) and 
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2(7-112), respectively. It follows from (100) that *GJ°(«) -* as 2-* 
oo, for n 2 <j<k 2 ,n 2 <l<k 2 . Hence, from (98), we have G,(z) -» as 
2 ^ oo, for n 2 <j<k 2 . It follows that we may expand Gj(z) in partial 
fractions in the form 

aw-Xerfe " ,a ' a * b (101) 

where ao = 1, and 

& r = lim [(2 - l/ffr)Gy(2)]. (102) 

Z—l/or 

From (93 , for »i < i < A:,, (94), and (101), we deduce that 

2./2+1 

Pu = - S <Tr +1 "'"&r, ni < i < A,, n 2 <> < A 2 . (103) 

r-0 

We proceed with the explicit calculation of gjr. Let 

rf = lim [(2 - lfo r )Gf\z)l (104) 

2—l/or 

Then, from (98), (102), and (104), we obtain 

*2 

Orgjr = X (niOrPn^l ~ dip n x -\,l) g jr 

+ [morPmfla - (ai + a 2 )Pn l s H -i]gp ) . (105) 

As we noted in Section IV, </>„ 2 (w r ) ^ 0. It is shown in Appendix C 
that 

Tj(l/o r )=^-, n 2 ^j^k 2 , (106) 

with wr given by (47). Hence, from (100), (104), and (106), since a = 1, 
we obtain 



(/) = _(^£| P'WWW _ w^ ( 10 7) 



„ , .V/"" 9 fr(0)fr(0) = _ to<0) 

*° W *n 2 (0)if(l) H(l) 

from (137). We note that H(l) ^ 0, since we have assumed that ai ^ 
Ai. We also obtain 



S. 



F 



(n 2 \ "' <fr/((0 r )<ft/(t0 r ) 

~^ 2 ) Ml-ff P )H'(ffrW«r)' 



P-(2| WWfWW , r=l, ...,2 92 + l, (108) 



where the prime denotes derivative. 

From (27), (103), (105), and (107), since a - 1, it follows that 

#(l)eo = 1 (nip niJ - aip ni -u) + n 2 p n ^ - (ai + a 2 )p n^-i. (109) 

l— n 2 
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If we substitute for />„,./ andp ni —u, n 2 < / < £2, andp ni ,n 2 -i from (17), 
and sum on /, and use (13) and (128), we obtain 

H(l)e = Yi dm[pm»»,-i(l + Pm, ai)s„ 2 (-p m , a 2 )^ q2 (p m ) 

m=0 

+ n 2 S ni {pm, ai)Sn. 2 (—pm, a 2 )(j)„ 2 {p m ) 

- (ai + a 2 )s ni (p m , ai)sn 2 -i(-pm, a 2 )</> n2 (pm)]. (110) 

After reduction with the help of (46) and some recurrence relations, 
we deduce from (110) that 

*2 

H{l)e = - 2 dm<f>n 2 (pm)[PmS ni (l + pm, ai)S^-l(l ~ p m , <Z 2 ) 

m=n 

+ (l\Sn x {pm, aJSnz-A-pm, d 2 )]. (Ill) 

But, summation of (37) from /" = to n 2 — 1 shows that the term on 
the right-hand side of (111) is zero, so that 

H(l)e = 0. (112) 

Since H(l) ¥* for a x ¥> Ai, it follows that e = 0, as was shown in 
Section IV. 

Next, from (27), (103), (105), and (108), we deduce that 

a r (l - O r )H'(Or)<t>n 2 (Ur)e r 

k 2 /n 2 V" 2 

= Y (a\Pn x -\J— n\O r Pn x .l)\ 1 <t>l(Ur) 

i-nt \a 2 J 

+ [(ai + a 2 )p„„„ 2 -i - n 2 a r Pn„n i ]<Sir h (iOr), (113) 

for r = 1, ••• , 2q 2 + 1. If we substitute for p„ x j and p nx -\j, n 2 < / 
< k 2 , andp n ,, n2 -i from (17), we obtain 

k 2 
o r (l — Or)H'(o r )e r = X d m [aiS ni -i(p m , a\) — /iia r s„,(p m , ai)] 

k 2 / \'-"i k 2 

-Sn 2 (-pm,a 2 ) X ( J */(pm)^/(«r)/^n,(Wr) + £ d m Sn x (p m , di) 

l-n 2 \a 2 J m=0 

•[(cti + a 2 )Sn 2 -i(-p m , a 2 ) - n-zOrS^-pm, a 2 )]^ ni (p m ), ( 114 ) 

for r = 1, • • • , 2^2 4- 1. The sums on / are given explicitly by (144), or 
by (145) if p m = w r . 

If we set a r = 1 in the right-hand side of (114), so that co r = from 
(47), it is found after reduction that the expression reduces to — H(l)e , 
as given by (111), and hence is zero, from (112). It follows that 
(1 — a r ) may be factored from both sides of (114). Consequently, we 
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may consider the limiting case Oi -► Ai, and obtain a finite expression 
for e 92 +i, even though a q2 +i -* 1 88 ffi -» A\. 

Now that we have determined the constants {e r ) in terms of the 
constants {d m }, it remains to determine the latter. These are deter- 
mined from (37) for < j < n 2 - 1, (39) for n 2 + 1 < /' < k 2 and (49). 
The redundancy of (39) for j = n 2 was discussed in Section IV, and 
follows from the fact that H(l)e = 0. Thus, the constants {d m } are 
determined from the boundary conditions at i = m, <y* < n 2 - 1 and 
at i = k u n.2 + 1 < y* < /e 2 , and from the normalization condition, the 
boundary condition at i = k u j = n 2 being redundant. The constants 
{d m } have to be calculated numerically. 

VII. AN ALTERNATE DERIVATION 

In this final section, we will give another derivation of the expressions 
for Co and e r , r = 1, • • • , 2q 2 + 1, given in (109) and (113). This 
derivation is more direct, but somewhat obscure, and was not evident 
until the results had been established with the help of the generating 
function approach. 

We make use of the representation (27) for the probabilities pa, 
which holds for m < i < ki and n 2 < j ' < k 2 , and the boundary 
conditions at i -= Hi, n 2 <> < k 2 , as given by (35) and (36). In particular, 
from (27) we have 

2*2+1 

Pn 1 J= £ M>y(Wr), n 2 <j<k 2 . (115) 

r-0 

Next, from (27) and (36), with the help of (47) and (134), we obtain 

292+1 p 

*4-W- I -fc(«r), n 2 +l<y<* 2 . (116) 



r-0 Or 



Also, from (27) and (35), with the help of (47) and (140), we find that 

2^2+1 g r 

ai/>„,-i^ + (ai + ajpn^-i - <*i E — *»*M. (117) 

r-0 Or 

From (115) to (117), it follows that 

2 (OiP„ 1 -lJ-/llffmPn 11 >)[ — I $/(Wm) 

+ [(ai + a 2 )p ni ^j-i — »2ffii»Piii^t]^»t( w '") 

IT *[(— -»*) 2 (?] tft+toM 

r Zo \ffram / y-nj \a 2 / 



>-"2 



2«2+l 
r-0 



- »2«rS(«l«)^»,(Wr) 



(118) 
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Now, from (47), 

(<o r - w m ) = (<r m - a r )( — ni ). (119) 

\O r a m ) 

But, for a\ j6 Au we know that a r ¥* a m for r =£ m. Hence, if we set 
p = w m in (143), and make use of (119) and (142), we obtain 

«l I L I —I 9fWm)9A a r) 

OrOm J j-m \Cl2/ 

= n2</»n 2 (w m )(/>„ 2 (to r ), r ¥■ m, r, m = 0, • • • , 2q 2 + 1. (120) 

This implies that all the terms on the right-hand side of equation 
(118), except the ones corresponding to r = m, are zero. 
Next, from (13), (20), (25), and (133), it follows that 

(1 - a)H(a) = orftar-iM - mofctfc), (121) 

and 

^-[(l-o)H(o)] 
da 

= I ni \ J[a2^^-i(<o) - /i2a</>'^(to)] - n 2 ^(w). (122) 

Hence, from (145) we obtain 

fe-" 1 ) X feT W^)] 2 -*k)] 2 

= ^> r )H[(l-a)H(a)]} . r-0,...,2* 2 +l. (123) 
But, a = 1 and w = 0, and H{a r ) = 0, r - 1, • • • , 2g 2 + 1, so that 

(>-«]L-(s ffW , jn... f%1+ L < i24> 

The expressions in (109) and (113) for e and e r , r » 1, • • • , 2q 2 + 1, 
respectively, follow directly from (118), if we set m = and m = 
1, • • • , 2q 2 + 1, respectively, and make use of (120), (123), (124), and 
(137). 

APPENDIX A 

Properties of the Eigenfunctions 

We define s,(A, a) by the recurrence relation 

(a + i + X)si(\, a) 

= a(l - 6,o)s,-i(A, a) + (i + l)s l+ i(X, a); s (X, a) = 1, (125) 
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for i = 0, 1, • • • . Thus, s n (\, a) is a polynomial of degree n in both A 
and a, and it may be related to a Poisson-Charlier polynomial. 5,6 
However, we will give here the properties of s„(\, a) that we will need. 
An explicit formula is 1 

«^<* •&$£>>• (126) 

where 
(A) = 1, (A)r - MX + 1) • • • (A + r - 1), r = 1, 2, • • • . (127) 
It was also shown 1 that 

(i + l)s l+ i(A, a) = asi(\, a) + As,(A + 1, a), (128) 

and 

s,(A, a) = Si(X + 1, a) - (1 - S.ote-ifA + 1, a). (129) 

From (129), it follows that 

£ s,(A, a) = s n (X + 1, a), (130) 

i-O 

and, from (128) and (130), we deduce that 

£ (n - i)s«(A, a) = (1 - S n0 )s n -i(X + 2, a). (131) 

We now turn our attention to the Chebyshev polynomials of the 
second kind, 7 Ui(x). They may be denned by the recurrence relation 

2xU,(x) = U l+l (x) + U^iix); U-^x) = 0, U (x) = 1, (132) 
for / = 0, 1, • • • . From (12) and (132), it follows that 

(a 2 + n 2 - p)*i(p) = cuVi+iip) + rc 2 ¥,-i(p); 

*_ x (p) = 0, * (p) = 1. (133) 

From (13) and (133), we deduce that 

[a 2 (l - 8jk 2 ) + n 2 - p]<t>j(p) = a 2 <t>j-i{p) + n 2 (l - Sjkjh+iip), (134) 
forj<k 2 . Since 7 



Ui 
it follows that 



IH 



= i e r ~ i , d35) 

r-() 



*, (0 > = H /2 U(^) = i K (136) 

VW \2Va2nJ r-o \a 2 / 
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Hence, from (13), we have 

&(0) = f jjjj . (137) 

If we set p = w r in (134), and also use (134) as it stands, we obtain 

(W r - p)(j>j(p)4>j(Ur) = a>[<f)j l(p)4/(Cl>r) "" $Ap)<l>J -l(w r )] 

- n 2 (l - 5y* 2 )[^>(p)4»y+i((Or) - $/+i(p)$/(«r)]. (138) 
If we now multiply (138) by (n 2 /a 2 ) J ~"', and sum on/, we deduce that 

(U r -p) I —J <j>j(p)<t>jM 

/-"a V" 2 / 

- Olfc*-! (p)6*(*>r) - <t>n.Xp)<t>n, -l(Wr)]. (139) 

But, from (21) and (22), 

[a 2 + /l 2 (l ~ a r ) - U r ]4>n.XUr) = /l 2 <f>n..,+ l(<Jr) , (140) 

for r = 0, • • • , 2q> + 1. Also, from (134), since q> > 1, 

(a 2 + n 2 - uv)4>„,,(wr) = a2^„,-i(a) r ) + n 2 4v,+i(wr). ( 141 ) 

From (140) and (141), it follows that 

a 2 0„,-i(io r ) = rc 2 a,4>„,(io,), r = 0, • • • , 2q 2 + 1. (142) 

Hence, from (139) and (142), we have 



(0) r -p) X 1 — 1 *>(p)*/(Wr) 

>-/.., \a 2 / 

= <t>nj,u r )[a 2 <t>n.,-\(p) - n>Or(t>n.,{p)], r - 0, • • • , 2g 2 + 1. (143) 

If we set p = pm in (143) and make use of the fact that p„„ m = 
0, ■ ■ • , k>, are the roots of equation (15), we deduce that 
k 2 /„ \J--2 

(tO r - p„,).S ;l .,(-p„„ a 2 ) £ (— I 0>(Pm)^/(o)r) 

= $n 2 (pm)<}>n 7 (u r )[a2Sn. r \(-p m , 2 ) — n 2 OrS„.,(-pm, Of)]. (144) 

Also, if we divide (143) by (co r — p) and let p — * u r , and use (142) and 
L'Hospital's rule, we find that 



'"hi \ a2 J 



= 4>nJ,Ur)[n>O r (}>'n.ilOr) ~ a 2 <$>'n. -l(Wr)], T = 0, • • • , 2q > + 1, (145) 

where the prime denotes derivative. 
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APPENDIX B 

Eigenvalues Corresponding to the Primary Queueing Region 

We consider here the zeros of 



where 



H(o) - U, - ^W(w) + n*frM, (146) 



w- |^-n, )(l -«). (147) 



As may be seen from (12) and (13), ^ 92 (to) and $„,(<<>) are both 
polynomials in to of degree q 2 . It follows that a q2+1 H{a) is a polynomial 
in a of degree 2g 2 + 1, so that H(a) = has 2q 2 + 1 roots. Now 7 

Ul{coB9) .^»l»i. (148) 

sin 



«m = a2 + «2 — 2vWi2 cosl ), m = 1, • • • , q 2 , (149) 



We define 

mm 
,q7+lj' 

and we note that < «i < • • • <a q . i . From (12), it follows that 

% 2 (<x m )=0, m = l,--',q 2 . (150) 

Next, we consider 

<*>„» m % 2 ( U ) - % rl (u), (151) 

and we note from (137) that <f>„ 2 (0) > 0. If q 2 = 1, then 4>„ 2 (ai) = -1, 
and the unique zero fii of <f>„ 2 (to) satisfies < fii < «i. For ^ 2 ^ 2, we 
define 

f s = a 2 + Ai2 - 2\/a 2 7i 2 cos( — I, s =• 1, • • • , q 2 — 1, (152) 

W 

and we note that 

< cti < Si < • • • < a 92 -i < & 2 -i < «7 2 - ( 153 > 

From (12) and (148), it follows that 

*Vi(k)-0, s=l, ..., 9 2-l. (154) 

Then, from (150), (151), (153), and (154), we deduce that 

4>„ 2 (ai) = -*« 2 -i(«i) < 0, <KtfO = **(fi) < 0, 
<M«2) = -¥«-i(a») > 0, • • • . (155) 

We let 

fe(Ai) =0, m - 1, .... ^2, (156) 
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with /?i < /? 2 < • • • <&q r Then, it follows that 

< /3, < a, < f i < /? 2 < a 2 < • • • < $Vi < /?,, 2 < a 9z . (157) 

Now, from (147), it is seen that to decreases from +» to as o 
increases from 0+ to min(l, a.\/ri\), and to increases from to +oo as a 
increases from max(l, ai/fli) to +oo. Hence y m and 8 m , m = 1, • • • , q 2 , 
are determined uniquely by the relations 

— - n\ 1(1 - y m ) = a m , 0< y m < mini 1, — J, 

j- - n, j(l - 8 m ) =$ m , 0<S m < minf 1, ^J. (158) 

It follows from (157) that 

< y q2 < 8g 2 < • • • < yi < 8\ < min(l, cti//ii). (159) 

Now, from (146) and (147), since <>„. 2 (0) > 0, we have H{a,\/n\) > 0. 
But, from (150) and (156) to (158), 

H(8 X ) = f #i, - j)% 2 (fii) < 0, H(yi) = n 2 4>n 2 M < °- < 16 °) 

Also, for q 2 > 2, 

H(« 2 ) = U, - flW <#») > 0, H(y 2 ) = n^fa) > 0, (161) 

and so on. Finally, consideration of the sign of ^^(to) for to —* +oo 
leads to 

ff<° + )=(:°°' ^ 92is T, n - d62) 

[ + oo, if q 2 is odd. 

From the above results, it follows that there are roots of H(a) = in 
the intervals (0, y„ 2 ) t (8 q .,, y ?2 -i), ••• , (62, Yi) an d (<5i, ai//ii). Since 
5i < 1, there are at least q 2 roots in the interval (0, 1). 
Next, from (159), we have 

(, a.\\ «i Ci a\ a\ ,.,„. 

1,— <— -< < ... <— — < . (163) 

n\) riiSi niYi ttiOgj n\y q . 

Also, from (146), (147), (150), and (156) to (158), we obtain 

H(-^-)=n i (l-8 i )%.,(p i )>0, tf(— )=/i 2 <Mai)<0, (164) 

and, for q 2 > 2, 

tf(-^r ) -re,(l- 8>)%,((1,)<0, H[— ) = n 2 <M« 2 ) >0, (165) 
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and so on. Hence there are q 2 roots of H(a) = in the intervals 



These q 2 roots are greater than 1, and we have seen that at least q 2 of 
the remaining 9-2 + 1 roots lie in the interval (0, 1). It is evident that 
the 2q 2 + 1 roots of H(a) = are positive and distinct. 

It remains to consider the root that lies in the interval (S u ai//ii). 
But, from (146) and (147), with the help of (13), we obtain 

ff(l) - %,(0)(Ai - «i), (166) 

where A 1 is as defined in (26). Also, we have shown that H(8i) < and 
H(ai/rii) > 0. Hence, if a t < A\ this root lies in the interval (Si, 1), and 
there are q 2 + 1 roots in the interval (0, 1), and q 2 roots greater than 
1. If ai > A\, then ai/ni > 1 and the root lies in the interval (1, fli/lti), 
and there are q 2 roots in the interval (0, 1), and q 2 + 1 roots greater 
than 1. If a\ = A h then one root is unity, and there are q 2 roots in the 
interval (0, 1), and q 2 roots greater than 1. 

APPENDIX C 

Results Pertaining to the Generating Function 

We derive here the solutions of (97), and we first consider I = n 2 . 
Then, 

[a 2 (l - 8 jk2 ) + n 2 - flGj^M 

= a*Gpl(z) + Ml - 8 Jk2 )G$l(z), (167) 

for n 2 + 1 </' < k 2 . It follows from (134) that 

Gj^(z) = G£ B> (s)*(f). n 2 <j<k 2 , (168) 

since 4»* 2 ({) = 1. If we sety = /i2 in (97), and substitute from (168), we 
obtain 



a 2 + n 2 [ 1 - - ) - f 



<M?) " «2^ 2+ i(nJGi n 2 2) (2) = - i. (169) 



With the help of (13), (95), and (133), we deduce that 



a 2 + « 2 ( 1 J - C 

- (l - i)[(ni - ais)« ft (t) + n 2 <M{)] = f 1 - jVd/«). ( 17 °) 

from (20) and (25). Hence, 

(1 - 2)rY(l/2)Gi" 2) (2) = 1. (171) 
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Since T„ 2 (z) = 1, from (99), we have established (100) for I - n 2 . 
Next, we consider / = k 2 . Then, from (97) 

a2 + n J 1 -^i\-^ 



(<*a) 



Gp\z) 



,(k 2 ) 



<*->, 



= a 2 (l - 8 jn2 )G y ;>!(z) + n 2 G%l(z), (172) 

for n 2 <j<k 2 - 1. But, it may be verified from (95), (99), and (133), 
that 



aa + nofl-aj-f 



T,(z) 



- a 2 {\ - SjnJTj-Az) + /i 2 T, +1 (2), (173) 

for./ > n 2 . It follows that 

G^iz) = G { n f(z)Tj(z), n 2 <> < k 2> (174) 

since T„ 2 (z) = 1. If we sety = k 2 in (97), and substitute from (174), we 
obtain 



[(n 2 - £)T* 2 (2) - a 2 T A2 - 1 (2)]Gi* 2) ( 2 ) = --. 



(175) 



With the help of (13), (95), (99), and (133), this may be written in the 
form 



ai 



- ) (1 - z)[(n x - a i2 )* q2 (n + n 2 <f> n . 2 (£)]G { n f(z) = 1. (176) 



It follows from (170) that 



a 2 



n 2 



- (1 - z)H(l/z)G { n ^(z) = 1. 



(177) 



Since <f>k 2 (£) = 1, we have established (100) for / = k 2 . 
Finally, we consider n 2 < I < k 2 . Then, from (97), 



«, + no(i-a| - c 



Gi'Hz) 



= a 2 {\ - 5 y „ 2 )GJ / J 1 (2) + n 2 G ( j'Uz), (178) 
for n 2 <j<l—l, and 

[a 2 (l - 8 jk2 ) + n 2 - nG} l) (z) = a 2 G?2 x (z) + n 2 (l - Sj k2 )G}'^(z), (179) 
for / + 1 <> < k 2 . It follows from (134) and (173) that 



G in (z) = {Gli:(z)Tj(z), n 2 <y</, 



(180) 
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since T„.,(z) = 1 and fa.XO = 1. The consistency of the representations 
fory = / requires that 

Also, if we set 7 = / in (97), and substitute from (180), we obtain 
(oa + n 2 - OG { "Az)<t»(0 ~ a 2 G%(z)Ti-i(z) 

-mGff<«)* + itf)--j. d82) 

It remains to solve (181) and (182) for G\ l ,l{z) and G^'tz). 
Now, since / < k>, it follows from (134) and (182) that 

1 



*[Gff«*-i<fl - G£U)T,-i(*>] - - -. 

If we eliminate Q%(z) with the help of (181), we find that 
a 2 2GJ£(z)[T/(e)<fc-,tt) - ^(flTi-ite)] = -T/(z). 
But, from (134) and (173), for n><j < k 2 , 

(a> + n> - $)<frA0 = Cfeto-itf) + Mtf+itf). 
and 

(a, + n> - f)Ty(z) = 02Ty-,(2) + n 2 T y+I ( 2 ). 

Hence, 

flk[Ty(s)#-i(f) - <M£>T/-l(*)] 

= n 2 [T, +1 (z)<M£) - 4g+i(flTj(fi)l (187) 
for n 2 <j< k 2 . Since n 2 < / < fe, it follows that 
Ti(*)*r-i(f) - */tf)T,-,<2) 



(183) 

(184) 
(185) 
(186) 



a. 



l-n>- 1 



- I - 1 PViW*4<fl - *%+i(f)T % («)]. (188) 



But, from (99) and (133), 
n 2 [T n2+ ,(z)^ 2 (f ) - ^, + .(f)T„ 2 (2)] 



a 2 + nil 1 I - ? 



<M?> - nrfn.+itf). (189) 



Hence, from (170), (188), and (189), we obtain 



«2|T,(2K„,(i)-<v(s)T, ,U)| (^) (l-iWij. (190) 
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It follows from (184) that 

^J (1 - z)h(^JG 1 1!(z) = T/( 2 ). (191) 

In view of (180), (181), and (191), we have established (100) for 
n.2< I < k 2 . 

We will now establish the relationship (106). From (21) and (22), it 
follows that 

[a 2 + Al 2 (l - Or) - Ur]^j(«r) = n 2 <f)n 2 +\((*)r). (192) 

Hence, from (47), (95), and (189), since T„,(z) = 1, we have 

T„ 2+ ,(l/a r )^(« r ) = (J>*, + i(a> r ). (193) 

As we noted in Section IV, <t>„.,(w r ) ¥= 0. Since T„., (l/o r ) = 1, we see that 
(106) holds for bothy - n 2 andy = n> + 1. But, from (47), (95), (185), 
and (186), <t>j((ti r ) and T/(l/a r ) are both solutions of 

{a 2 + n 2 — u r )hj = a 2 hj-i + n 2 h j+i , n 2 <j < k 2 . (194) 

It follows that (106) holds for n 2 <y < k 2 . 
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